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The free energies of the FCC, BCC, HCP and Simple Cubic phases for hard spheres are calculated 
as a function of density using the Fundamental Measure Theory models of Rosenfeld et al (PRE 



^ ; 55, 4245 (1997)), Tarazona (PRL 84, 694 (2001)) and Roth et al (J. Phys.: Cond. Matt. 14, 

^ ■ 12063 (2002)) in the Gaussian approximation. For the FCC phase, the present work confirms the 

"O - 

p^ . vanishing of the Lindemann parameter (i.e. vanishing of the width of the Gaussians) near close 

o, 

O , packing for all three models and the results for the HCP phase are nearly identical. For the BCC 
phase and for packing fractions above rj ~ 0.56, all three theories show multiple solid structures 

differing in the widths of the Gaussians. In all three cases, one of these structures shows the 

00 . 

Tj- , expected vanishing of the Lindemann parameter at close packing, but this physical structure is 

^D ' only thermodynamically favored over the unphysical structures in the Tarazona theory and even 

^ '. 

^P ■ then, some unphysical behavior persists at lower densities. The simple cubic phase is stabilized in 



the model of Rosenfeld et al. for a range of densities and in the Tarazona model only very near 
close-packing. 

PACS numbers: 64.10.+h, 05.20.-y, 05.70.Ce 



I. INTRODUCTION 

The Fundamental Measure Theory (FMT) approach to building approximate free energy 
density functionals has proven very successful in describingthe properties of inhomogeneous 
hard-sphere fluids including the hard-sphere FCC solid[l|, |2, Isl y, IS Isf . The original form of 
FMT proposed by Rosenfeld[l|,l4| was based on a generalization of the ideas underlying scaled 
particle theory. The Rosenfeld functional gave a useful description of an inhomogeneous fluid 
but did not stabilize the FCC solid. One important formal property of the functional was 
that its second functional derivative with respect to the density reduced to the Percus- 
Yevik direct correlation function in the uniform liquid limit. Subsequent refinements were 
motivated in particular by the requirement that the free energy functional reduce to the 



Q. 



known exact form when the density was restricted to zero- and one-dimensional systems 
The resulting theories, discussed in more detail below, retains all of the useful properties 
of the original Rosenfeld functional while also stabilizing the FCC solid and predicting the 
liquid-solid coexistence reasonably accurately [sl, 14 y, |y|- The theory also gives a good 
description of the mean-squared displacement of the atoms in the solid phase including 
the nontrivial prediction that the mean-squared displacement goes to zero as the density 
approaches close packing. 

The purpose of this paper is to present the results of the application of the FMT to other 
crystalline structures of the hard-sphere solid. There are two good reasons for carrying out 
such a study. First, from a practical point of view, one would like to be able to use the FMT 



to study solid-solid phase transitions as well as a basis for the study, via t 



lermo dynamic 



perturbation theory, of non-FCC solids along the lines discussed in refs. JT], |8|, |9|, |lO| . To do 
so with any confidence first requires that the predictions of the FMT for non-FCC hard- 
sphere solids be well understood. The second reason is that more loosely-packed crystal 
structures may provide a more demanding test of the theory than does the FCC phase since 
the structure of the FCC phase is not that dissimilar to that of the liquid: other structures 
differ more from the liquid structure, e.g. in terms of the nearest-neighbor coordination. As 
shown below, this expectation is borne out and the predictions of the properties of the BCC 
phase in particular are not as satisfactory as for the FCC phase. 

The remainder of this paper is organized as follows. Section |n] reviews the elements of 
Density Functional Theory (DFT) and FMT and discusses the difference between the three 



FMTs considered in this work. The results of the calculations are presented in Section IIIII 
where it is shown that the theories give very different results for different lattice structures. 
The final Section summarizes the results of these calculations and discusses obstacles that 
prevent any simple modification of these theories to give a better description of the BCC 
phase. 

II. FUNDAMENTAL MEASURE THEORY 

As is usual in DFT, the Helmholtz free energy F is written as a sum of an ideal contri- 
bution Fid and an excess contribution, Fex- The former is given by 

(5F,d [p] = I [p (T") In A'p (r ) - p (r )] dr (1) 

where p{r ) is the local density, N = Jp{r)dr is the total number of particles and 
P = l/ksT is the inverse temperature. In the FMT approximation, the excess free energy 



functional is written in terms of a set of local functions 

pf,4p] = /"^/50,(K(r,[p])})rfr (2) 

which, as indicated, include a functional dependence on the density through a family of local 
functionals of the form 



ria {~r, [p]) = j Wa{~r - ~r') p {~r') d~r'. (3) 

The weights w^ ("r*) occurring in eq.© are 6 (| — r), 5 (r — |) , fS (r — |) and rr5 (r — |) 
yielding the local density functionals which will be denoted ri{r),s{r), v {r ), T {r ) 
respectively. For a uniform liquid, in which p{r ) = p, one has, in three dimensions, that 
1] {r ) = ^na^J) which is the usual definition of the packing fraction. The names of the other 
quantities are motivated by their scalar, vector and tensor natures respectively. The first 
two contributions to the excess free energy are 

s (T^) \n {1 - 1] (T^)) (4) 



TTCr^ 



2-na {1 — 1] (^)) 
These are the same as in the original Rosenfeld theory. The third function has been the 

focus of most efforts to refine the FMT and we write it in the form 

Sn [1 — r] [ r )) 



10. 



where the only fundamental constraint on / is that linip^o / = 1 H • Three proposals in 
the literature, aside from the original form / = |s^ {^ ), which does not stabilize the bulk 
FCC solid, are 

fnsLT = lsHr){l-eir))\e = p^^ (6) 

fE = l{i-v {^)f (/?/L iv (^)) + log (1 - V (^)) - Y^^) fT 



where Jrslt was proposed by Rosenfeld et al|33 /t is the tensor form proposed by Tarazona 
and Je is a heuristic modification of fx that allows for the insertion of an empirical equation 



of state for the liquid, Pf^xiv)- I^ t 
be the Carnahan-Starling equation 



:;he latter case, if the equation of state is chosen to 
ll|, the resulting free energy functional is the single- 
species version of the so-called "White Bear" functional of Roth et al.[6], which was also 
proposed by Tarazonapj. The first and second forms are closely related. In fact, writing 
the tensor density as T (T^) = U (T^) + |s (T^) 1 , where U (T*) is traceless, a natural 
approximation for the traceless part is to make it proportional to the only traceless tensor 
that can be formed from the other densities 

U [r)^ A -^ (7) 

where A is a constant and the denominator is needed so that U {r ) scales linearly with the 
density in the sense that p (T^) -^ \p (T^) implies that U (T^) ^ \U i/r^)- Substituting this 
into fx gives fx = |s^ (T*) (1 — 3^^ + 3^4^'^ — A^C^^) which reduces to /rslt provided that 
one takes A = 1. The choice A = corresponds to the original FMT theory of Rosenfeld. 
Note that one could also divide by \v ( r ) | rather than s ( r ) in this ansatz, which is 
equivalent to the replacement A — > A'^~^ ("r^) giving fx = |s^ (T^) (1 — 3^^ + A' {3 — A''^)^^). 
Taking A' = 1 gives the "interpolation" form for 03 discussed in ref. J3|. 

The calculations reported here were performed using the standard approximation of the 
density as a sum of Gaussians 

n 

where / (r) = xq (f ) exp (— ar^). For a perfect lattice, one has occupancy xq = 1 whereas 
values Xq < 1 characterize a solid with vacancies. The parameter a characterizes the width 



of the Gaussians and < i?„ ^ is the set of lattice vectors of the desired crystal structure. The 
density can also be written in terms of components in reciprocal space as 



p[r 



■) = ^ exp [iKn ■ T'j pn (9) 



where < Kn \ is the corresponding set of reciprocal lattice vectors and p„ = 
xoNceiiPceiiGW {~Kn/^(^) with Nceii being the number of lattice sites per unit cell and 
Pceii being the number of unit cells per unit volume. Thus, NceiiPceii is just the num- 
ber density of the perfect solid. From the fundamental theorems of DFT, the equilib- 
rium values of the parameters xo,Pceii and a are those which minimize the grand potential, 
Pfi = [3Fid [p] + (3Fex [p] — (3pN for fixed chemical potential p and volume V . The number 
of atoms is A^ = JiV where the average density \s J) = ^ jy p{ r )d r = x^NceiiPceU- 

The spatial integrals in equations (0)-(j21) were performed by sampling the arguments on 
a grid of points covering the standard unit cell. The primary consideration in fixing the 
number of points in the grid is that the spacing be sufficiently fine so as to be able to sample 
the width of the Gaussians, (Ai?) = ^. In the calculations discussed below, the spacing 
of the points in each Cartesian direction was set to of dr = a/{2q) where a is the lattice 
parameter and q was initially set to 20 and then doubled until dr/AR < 0.5. Increasing the 
number of points beyond this limit had no significant effect on the results. Note that in the 
Gaussian approximation and in the limit of aa^ >> 1 , the ideal contribution to the free 
energy approaches the asymptotic value 

/J/«-x,(|l„(^)-2,5). (10) 

This approximation was used for aa"^ > 200. The plots of the free energy as a function 
of a presented below were all generated based on sampling the free energy at intervals of 
6aa^ = 1,5,10,50 and 250 for acr^ < 20,200,2000,10,000 and 20,000, respectively. The 
free energy minima were determined using the the Simplex algorithm of Nelder and Mead, 
see e.g. ref. |l2|, as implemented in the Gnu Scientific Library (l3l|. which was terminated 
when the simplex size was smaller than 10^^. 
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FIG. 1: Dimensionless Helmholtz free energy per atom, (if = f3F/N , as a function of a for the 
FCC hard-sphere sohd as calculated using (a) the RSLT theory, (b) the Tarazona theory and (c) 
the White Bear theory. From bottom to top, the curves are for r/ = 0.56 — 0.68 in steps of 0.02. 
The arrows indicate the positions of the secondary solid minima. 

III. RESULTS 

A. The FCC crystal 



Figure Q shows the free energies as a function of a for the FCC solid as calculated from 



the three theories for a variety of average packing fractions r] = ^pcr ■ These calculations 



are performed for Xq = 1 which is known to be correct at the solid niininiuni|3| and which 
we have verified is very close to the free energy minimum at all values of a. All three 
theories show a minimum for all densities at a = corresponding to the liquid phase. This 
is unphysical for the higher densities which are well above random close packing, fj^^p ~ 0.64, 
and is a result of the use of the Percus-Yevik and Carnahan-Starling approximate equations 
of state which are only singular at r/ = 1. On the other hand, note the consistency in the 
location of the free energy minimum at finite a as well as of the value of the free energy at 
the minimum. While the location of the minimum in RSLT theory is somewhat higher than 



for the other two theories, the differences are not large. For all three theories, the location 
of the solid minimum increases rapidly with density which is a characteristic feature of the 
FMTJ3|,|5|. The Tarazona and White Bear functionals differ primarily in the small-a region, 
due to the difference in the equations of state of the liquid, but at large a the differences are 
negligible. The main difference between the tensor theories and the RSLT theory is also in 
the small-a region where, for r/ > 0.64 the RSLT theory predicts much larger energy barriers 
between the liquid and solid minima than do the tensor theories. Note also that for a narrow 
range of densities, 0.61 ^f]^ 0.63, the RSLT theory actually shows two minima at non-zero 
values of a. This appears to be quite unphysical since all atoms are, by hypothesis, assumed 
to be identical and so it is hard to see how there could be two metastable states with the 
same density but different vibrational amplitudes. 

The equilibrium free energy of the solid phase as a function of average density was de- 
termined by minimizing the free energy functional with respect to both a and the lattice 
parameter. From this, the densities of the coexisting liquid and solid as predicted by each 
theory were determined in the usual way by finding liquid and solid states with the same 
pressure and chemical potential. The RSLT theory gives the coexisting liquid and solid at 
packing fractions 0.491 and 0.534 respectively (compared to 0.491 and 0.540 reported in |3| 
and 0.491 and 0.534 reported in 14 1), the Tarazona theory gives 0.466 and 0.516 (compared 
to 0.467 and 0.516 reported in ref. 5 ) and the White Bear theory 0.489 and 0.535 (compared 
to 0.489 and 0.536 reported in ref. 5|). For comparison, the values from simulation are 0.492 
and 0.545 [l^. The reason that the Tarazona theory gives relatively poor predictions for 
coexistence is, paradoxically, because it gives a rather good description of the solid while 
the liquid equation of state is still that of the Percus-Yevik theory which is inaccurate at 
the high densities that characterize coexistence '5|. This is one of the main motivations for 
introducing the empirical equation of state so that a free energy functional is obtained that 
yields accurate values over the entire range of structures. 

The predicted Lindemann parameter, L which is the ratio of the root-mean-square dis- 
placement to the nearest neighbor distance, is shown in FigI21 For the FCC crystal, it 
decreases monotonically with increasing density and extrapolates to zero at a density near 
the FCC close-packing density of TipQ(j_(jp = 7r-\/2/6 ~ 0.74. 

It is worth noting that similar calulations have been performed for the closely related 
HCP phases. In all cases, the results closely mirrored those for the FCC phase shown here. 
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FIG. 2: The Lindemann parameter for the FCC phase as a function of rj as calculated using the 
RSLT theory (lower line) , the Tarazona theory (middle curve labeled T) and the White Bear 
theory (upper curve labeled WB). Also shown as dotted lines are the quadratic interpolation of 
the curves to L = based on the data for rj > 0.60. 

In particular, the free energies were always larger than, but within 1% of, the FCC free 
energy. 

B. The BCC crystal 



Figure El shows the same calculations for the BCC solid. In this case, all three theories give 
multiple solid minima over some range of densities. For f] > 0.56,the RSLT theory gives two 
solid minima with the low-density minimum having lower free energy than the high-density 
minimum and thus being the preferred state. The Tarazona theory shows multiple minima 
for f] > 0.61 with energies that differ by less than one percent. In the range 0.61 < rJ < 0.63 
the low-a branch has slightly lower free energy while for 0.635 < fj the high-a minimum is 
favored. The White Bear functional only develops a high-a minimum for 0.63 < rj and in 
all cases the low-a minimum has lower free energy. 

As seen in Fig. El the behavior of the Lindemann parameter for the BCC crystal is more 
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FIG. 3: Same as ^ for the BCC phase. In (b) and (c), the arrows show the 
positions of the sohd minima. From bottom to top, the curves are for r\ = 

0.52, 0.54, 0.56, 0.58, 0.59, 0.60, 0.61, 0.62, 0.63, 0.64, 0.65. 



complex than in the case of the FCC phase. All three theories give a "low-a" branch which is 
not monotonic: the Lindemann parameter decreases with increasing density at low densities, 
but then increases with increasing density at high densities. This behavior is similar to that 
of the effective-liquid-type DFTsj^, ll^ and is considered quite unphysical since, if the trend 
persisted up to close packing, it would imply that spheres are able to interpenetrate. Rather, 
one expects the mean-squared displacement to decrease uniformly with increasing density 
as in the case of the FCC crystal. All three FMTs also give a second, "high-a" branch 
which does behave physically at high densities and, in particular, the Lindemann parameter 
extrapolates to zero near the BCC close packing density of fJBcc-CP — ^rv^/S ~ 0.68. 
However, for with the RSLT theory, up to r^ = 0.625, the well-behaved branch has higher 
free energy than does the unphysical branch. The Tarazona theory gives similar results 
except that the high-a minima are thermodynamically favored at high density so that, 
as noted by Tarazonaj^ |5|, the prediction of the theory is that the Lindemann parameter 
appears to vanish near close-packing, as it should. Nevertheless, the prediction of the theory 
is that the Lindemann parameter increases with increasing density at intermediate densities 
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FIG. 4: The Lindemann parameter for the BCC phase as a function of packing fraction, r], as 
calculated using the RSLT theory, the Tarazona theory (labeled T) and the White Bear theory 
(labeled WB). Both the low-a and high-a branches are shown with the stable branch being drawn 
with full lines and the unstable branch with dashed lines. Also shown as dotted lines are the 
quadratic interpolation of the curves to L = based on the data for rj > 0.60. 

as well as giving a discontinuous jump in the Lindemann parameter for a packing fraction 
near 0.60, all of which seem somewhat unphysical. Finally, the White Bear theory also has 
a second branch but it only appears at very high densities and also has higher free energy 
than does the low-a branch. Comparison of Figures 2b and 2c shows that the reason for 
the observed behavior is that the Carnahan-Starling equation of state used in the White 
Bear functional lowers the liquid free energy relative to the Percus-Yevik equation of state 
that comes out of the Tarazona theory. As a increases, the effect lessens until at high a 
the theories are almost identical. Thus, the free energy of the low-a minimum is somewhat 
lower in the White Bear theory leading to its stability persisting even up to high densities. 
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FIG. 5: Same as Fig. ^for the SC phase. From bottom to top, the curves are for r/ in the range 
[0.43,0.52] in steps of 0.01. 

C. The SC phase 

For a purely repulsive potential, the simple cubic phase is expected to be unstable with 
respect to shear. Nevertheless, there is no reason that DFT should not give a well-defined 
free energy for this structure since it should be metastable. Indeed, at least one of the 
effective-hquid theories, the GELA, does give a stable structure with small values of a|lfi|. 

Figure |31 shows the free energy curves calculated using the three theories. The results 
from the RSLT theory are much as expected: the simple cubic phase stabilizes at quite low 
densities, 1] ~ 0.45, and the location of the free energy minimum increases with increasing 
density and appears to diverge as the simple cubic close packing, fjgQ_(jp = n/6 ~ 0.524 is 
approached. In this case, there is no indication of multiple minima. In contrast, the tensor 
theories show no sign of a a 7^ minimum except at densities very close to close packing. 
Extending the calculations up to aa'^ = 20, 000, no minimum was found in the rj = 0.51 
curve using either of the tensor models while for fj = 0.52 the Tarazona theory stablizes the 
solid at aa"^ ~ 11, 500 and there is again no minimum using the White Bear theory. 
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TABLE I: The minimum packing fractions at which the various FMT models predict the different 
crystal structures become stable and the approximate value of the Gaussian parameter a at this 
density. 



Structure 



Theory 



'^sol 



ao 



FCC 



BCC 



SC 

sc 

SC 



RSLT 

T 

WB 

RSLT 

T 

WB 

RSLT 

T 
WB 



0.49 
0.46 
0.48 
0.51 
0.48 
0.51 
0.47 
0.52 



27 
22 
26 
14 
20 
24 
1150 
11,250 



D. Relative Stability 

The last question we adress is the relative stability of the various structures. Table E] 
shows the minimum densities at which each structure becomes stable in the various FMT 
models. For the FCC and BCC structures, the RSLT and WB theories give quite similar 
results while the Tarazona theory generally stabilizes the solid at a lower density about 5% 
below the other theories. In all cases, the FCC phase is stable at lower densities than the 
BCC phase. 

The free energies as functions of the average density are shown in Fig. IHl In all models, 
the FCC phase has the lowest free energy of any solid phase at all densities. The free energy 
of the BCC phase is very close to that of the FCC at low densities and steadily diverges 
at increasing densities. The FCC phase has slightly higher free energy than the liquid at 
the lowest densities. The low-a branch of the BCC phase always has lower free energy than 
the liquid, although in the RSLT theory they are almost identical. In the RSLT, the simple 
cubic structure has much higher free energy than any other phase. 
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FIG. 6: The free energy minima for the hquid phase and the various crystal structures as functions 
of the average packing fraction calculated using (a) the RSLT theory, (b) the Tarazona theory and 
(c) the White Bear theory. The lowest curves are for the FCC structure, the intermediate curves 
are the low-a (full) and high-a (dashed) BCC branches and the light curve that spans the width 
of the figures are for the liquid phase. In (a), the upper curve is for the simple cubic structure. 

IV. CONCLUSIONS 

In summary, the present calculations confirm the fact that the three FMT's discussed all 
give reasonable results for the FCC phase at all densities up to FCC close-packing with the 
only exception being the appearance of a spurious low-a minimum at high densities in the 
RSLT theory. On the other hand, when applied to the BCC phase, these results indicate 
several unphysical properties of the Fundamental Measure Theory functionals which must 
be carefully considered when using them, e.g., as the basis for thermodynamic perturbation 
theory. All of the theories appear to give reasonable results at moderate densities, say 
up to r/ ~ 0.56. However, the simplest model, the RSLT, gives multiple solid minima 
for all packing fractions in the BCC phase above about 0.56. The Tensor theories do not 
give multiple minima for the FCC phase but they do for the BCC phase. As the density 
increases, the stable minima in the Tarazona theory switches from the unphysical branch 
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to the physically well-behaved branch whereas the White Bear theory always picks out the 
unphysical branch. Only the RSLT and Tarazona theories stablize the simple cubic phase, at 
least for aa'^ < 20, 000, but the Lindemann parameter is, in this case, well-behaved and there 
are no spurious minima. Unfortunately for the goal of choosing a single theory to describe 
both the liquid and all solid phases, the RSLT and WB theories give poor results for the 
BCC phase Lindemann parameter while the Tarazona theory gives the worst quantitative 
description of FCC-liquid coexistence. 

It is interesting to speculate on the resolution of the common deficiencies in these the- 
ories. It seems reasonable to attribute the unphysical behavior at high densities, which is 
manifested by the fact that the unphysical low-a minimum is favored, to a poor description 
of the liquid phase: the free energy of the liquid must clearly diverge at high densities and 
in fact it has long been debated whether the hard-sphere liquid-phase equation of state pos- 
sesses a singularity at high densities 17|] - at, e.g., random close packing[l8| or FCC close 
packing |l9l|. If the free energy of the liquid did indeed diverge at random close packing, 
then the curves shown in Figs. Q and El would be deformed with the low-a free energies 
"pulled up" . This could conceivable eliminate altogether the low-a minima and thus give a 
reasonable picture of the BCC phase. However, as appealing as this scenario is, it cannot be 
implemented naively using eq. © to give the equivalent of the "White Bear" improvement to 
the Tarazona theory. This is because eq. (0) requires that the equation of state be evaluated 
at the local packing fraction r/( r ) and for the solid phase there are always points at which 
?7( r ) is very nearly equal to one unless a is very small. (In particular, it is easy to see that 
liniQ^oo ''?( ) = 1.) Thus, simply inserting an equation of state with a divergence for fj < 1 
is not an option since this will lead to divergent free energies in the solid phase, even for the 
already well-described FCC phase. 

One approximation not tested here is the use of the Gaussian profile for the densities. 
Tarazona has performed calculations using more general parameterizations of the density 
and he reports that the generalization beyond the Gaussian is of negligble imortance for 
the equation of state[4]. Groh and Mulder have also supplied evidence for the FCC phase 
that the Gaussian approximation is quite accurate, especially near the peaks of the density 
distribution l20l|. This suggests that the effects reported here are not due to the Gaussian 
approximation but definitive proof will require further calculations using non-Gaussian pro- 
files. 



14 



In summary, the Fundamental Measure Theories work well for the FCC phase at all den- 
sities while for the BCC phase, the FMTs are best applied to the solid phase for densities 
near liquid-solid coexistence. While they can in principle be used to model the BCC phase 
at higher densities, and while they do show some desirable features at high density like the 
vanishing of the Lindemann parameter at close packing, they also suffer from the same prob- 
lem that plagued earlier DFT models: namely, the unphysical behavior of the Lindemann 
parameter for the BCC phase at intermediate densities. On the other hand, the SC phase 
is well modeled, at least using the RSLT theory, perhaps because it only exists at relatively 
low densities. 
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